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Abstract. This article discusses the classical problem of how to calculate r n (m), the number of ways to 
represent an integer to by a sum of n squares from a computational efficiency viewpoint. Although this 
problem has been studied in great detail, there are very few formulas given for the purpose of computing 
r n (m) quickly. More precisely, for fixed n, we want a formula for r n (m) that computes in log-polynomial 
time (with respect to m) when the prime factorization of m is given. Restricting to even n, we can view 
n (q), the theta function associated to sums of n squares, as a modular form of weight n/2 on Ti(4). In 
particular, we show that for only a small finite list of n can 6 n be written as a linear combination consisting 
entirely of Eisenstein series and cusp forms with complex multiplication. These are the only n that give 
rise to "elementary" formulas for r n (m), i.e. formulas such that for a prime p, r n (p) can be calculated 
in C(log(p))-time. Viewing 6 n (q) as one of the simpler examples of modular forms that are not strictly 
Eisenstein, this result motivates the necessity of a log-polynomial time algorithm that directly calculates the 
Fourier coefficients of modular forms in the generic situation when there is no such formula, as described in 
Couveignes and Edixhoven's forthcoming book (for level 1 cases) and Peter Bruin's Ph.D. thesis (for higher 
level, including 4). 

Introduction 

We are interested in the even case of the classical "sums of squares" problem. For positive n € Z and positive 
to e Z, let 

r n (m) := #{x = (xi, ...X n ) € Z" : x\ + x% + ... + x n = m}. (1) 

We wish to understand r„(m) by analyzing the generating function attached to r n (m), denoted by 9 n . 
As a complex function on the open unit disk T> C C, 9 n is the theta function attached to the standard 
rt-dimcnsional quadratic form, 



e n (q)= f> n (TO)<f" = l + 2n. g + 4Q V+sQ ■ q 3 



. I " I +2// 



rn—0 

Equivalently, one can define n (q) by using the multiplicative property of r n (m 

0n(q)=0i(q) n = (l + 2q + 2q 4 + 2q 9 + ...) n . 



■q 4 + .... (2) 
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Furthermore, if we view T> as the image of the upper half plane i = {z e C : Im(z) > 0} under the map 
z i— > e 27 " z , then we may write, 

oo 

9 n (z) = J2 r n (m)e 2 * lmz . 

It satisfies the equations (see [IT], 3.2), 

0„(-l/4z) = (2z/z)"/ 2 n (z), 0„(z + 1) = e n (z). (3) 

Since we restrict ourselves to the case where n is even, there is no need to choose a square root. The above 
symmetric properties illustrate that 9 n (as a function of the H-coordinate z) is a modular form of weight 
k = 5 on the congruence subgroup Ti(4), consisting of matrices 7 € SL2(Z) such that 7 = (J J ) (mod 4) 
(see |B], 1.2 or [T7], 3.2). 

We wish to analyze for which n does #„(z) have coefficients that are "straightforward to compute", i.e. 
for fixed n, when r n (m) can be written as a function on m 6 Z such that each Fourier coefficient a m = r n (m) 
of 6 n {z) can be computed in 0(log(m))-time, given the prime factorization of m. Equivalently, for a prime 
p, we want formulas for r n (p) that can be computed in time polynomial to the number of digits of p. We 
therefore establish the following definition. 

Definition 1. A modular form / of integral weight is elementary if and only if / is a linear combination of 
Eisenstein series and cusp forms with complex multiplication (CM cusp forms) as defined in Section [Tj 

Denote the space of modular forms on a congruence subgroup T of weight k as Mk(T), and let its subspace 
of cusp forms be <Sfc(r). The orthogonal complement (with respect to the Petersson inner product) is denoted 
£fc(r), the space of Eisenstein series. Finally, we define the subspace S^ m (T) C <Sfc(r) as the CM subspace, 
i.e. the space generated by those cusp forms which are invariant under twisting by a quadratic character. 

The definition of elementary comes up quite naturally in a variety of places. In particular, elementary mod- 
ular forms are in correspondence with 2-dimensional potentially abelian £-adic representations of Gal(Q/Q), 
i.e. after passing to an open subgroup (equivalently, to a finite extension of Q), the representation factors 
through an abelian quotient. These turn out to be particularly nice examples of geometric representations 
(see [8], §6 for details). 

Furthermore, in [24], Serre calls cusp forms lacunary if the arithmetic density of the nonzero coefficents 
in the g-expansion is zero, and proves that cusp form / g Sk(T) is lacunary if and only if / £ S^. m (T). While 
it is false that 6 n is lacunary for any n > 4 (due to Lagrange), it is true that 9 n is elementary if and only if 
the cuspidal part in its decomposition are lacunary, i.e. contribute to the value of r n (m) very rarely from a 
density viewpoint, making most coefficients of 9 n a sum of Eisenstein coefficients. The following theorem is 
our main result. 

Theorem 1. Suppose n is even. Then 9 n is elementary if and only if n = 2,4,6,8, or 10. 

For the listed exceptional cases, the following formulas for the Fourier coefficients r n (m) were produced 
originally by Jacobi (for n — 2,4,6,8) and Liouville (n = 10) in the 19th century (see [12] and [14]). Let Xi 
denote the nontrivial Dirichlet character of conductor 4 extended multiplicatively to all positive integers. 
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ra(m) = 4\J2xi(d) 



r 6 (m) 



, dim 











r±{m) = 8 


E d 


+ 16 


E d 




d odd 

^ dm / 




d odd 



16 (E^(7) d2 |- 4 (E^ 2 



(4) 



(5) 



(6) 



7-s (m) 



rio(m) 



16 (E^ 3 ~ 32 E d3 + 256 E d3 



i d\m 




E * 4 

aez[i] 
\Nm(0)=m / 



(7) 



(8) 



It has been known that for even n < 10, the theta function is purely Eisenstein, and therefore elementary 
(see [20]). For even n > 8, Liouville, Glaisher, and Ramanujan tried to give formulas in the same vein as those 
given for r2,r4,rg, and rg; however, they began to involve various arithmetic and analytic functions for the 
cuspidal coefficients which were exceedingly difficult to compute, even for small m (see [9] and [E]). General 
formulas for all even n were also remarkably difficult due to the contribution from the cusp form, though 
the study of the Eisenstein part of n (q) allowed for understanding of the general growth of these functions 
(see [16j). However, as we will show, CM cusp forms have coefficients given by a Hecke character on an 
associated imaginary quadratic field, hence it is still possible to have log-polynomial computable coefficients 
when the cuspidal part of 9 n (z) has CM, i.e. elementary formulas can exist for r n (m) when n > 8. To prove 
the main theorem, we compute the dimensions of S n /2(Ti(4)) @ S^J 2 (T . It turns out that the first cusp 
form on Ti(4) is of weight 5 and is invariant under twisting by the quadratic character attached to Q(z); in 
fact, 5| m (ri(4)) = £5(171(4)). Thus, we can already conclude that 9 2 , O4, 9$, 9s, and 9 W are elementary since 
they are modular forms on Ti(4) of weight less than or equal to 5. 

For the converse, we first show that 5^™ 2 (ri(4)) is one-dimensional when § = 1 mod 4 and trivial oth- 
erwise. Since the Riemann-Roch formula implies that £ n / 2 (ri(4)) has dimension at most 3, we only need 
to show that the modular form 9 n is not in the span of a basis of (at most 3) Eisenstein series and 1 CM 
eigenform (if ^ = 1 mod 4). This can be done by constructing bases for f„/ 2 (Fi(4)) ® 5^™ 2 (ri(4)) and com- 
paring coefficients. Thus, 9 n is only elementary in the exceptional cases where the entire cuspidal subspace 
of X 5 (ri(4)) is CM. Note that for n > 10, 5^ 2 (ri(4)) C 5„ /2 (r 1 (4)). 

We can also prove the backwards direction of Theorem[T]by looking at the action of certain Hecke operators 
on the CM cuspidal subspace. In particular, we can show that the Hecke action on the cuspidal part of 9 n (q) 
is not the same as it would be on an element of 6>^" l 2 (ri(4)) when n > 10. 

Section 1 describes cusp forms with complex multiplication and pertinent properties, including their 
relation to Hecke characters of imaginary quadratic fields. In Section 2, we show that 9 n (q) is elementary 
for even n < 10 by a dimension argument. In Section 3, we show that 9 n (q) for n > 10 cannot be elementary 
by comparing it with the Eisenstein and CM cuspidal subspaces of A^ n /2(r!(4)) concluding the first proof 
of Theorem [1] In Section 4, we produce another proof of the fact that 9 n (q) is not elementary for all n > 10 
and even by contrasting the action of the Hecke operator T 3 on the CM cuspidal space with its action on 
9 n . Finally, in Section 5, we discuss the context of this result and the notion of elementary in modern 
computations of Fourier coefficients of modular forms. 
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1 Cusp forms with complex multiplication 

We recall the theory of CM cusp forms, and describe their relation to Hecke characters of imaginary quadratic 
fields. For a more detailed explanation, see |21) . 

Fix positive integers k and N, and let e be a Dirichlet character on (Z/iVZ) x such that £(—1) = (— 1) . 
Recall that Sk(To(N),e) is the e-eigenspace of the diamond operator (d) where d £ (Z/iVZ) x (recall that 
r (JV)/ri(iV) S (Z/NZ) X and the diamond operators represent the action of T (N) on 5 fc (Li(A^))). For a 
normalized eigenform / £ Sk(To(N), e), let Kf be the number field generated by its Hecke eigenvalues. Note 
that Kf is either totally real or has complex multiplication, and it contains the image of the Nebentypus e. 
It is real if and only if e factors through {±1} C C x and 

e(p)a p — a p V primes p \ N, 

where a p denotes the coefficient of q p in the Fourier expansion of / (see Prop 3.2, [21]). 

For a newform / £ Sk(To(N),e), we can twist by a Dirichlet character ip defined on (Z/DZ) X as follows: 

oo 

f®<p = J2 f(n)a n q n £ S k (T (ND 2 ),ep 2 ). 

n=l 

Moreover, f <£) ip is an eigenform since the action of the Hecke operator T p for p\ ND has eigenvalue ip(p)a p . 

An eigenform / has complex multiplication (or CM) by ip if / <£> tp = f. One must check that ip(p)a p = a p 
(equivalently, either ip(p) = 1 or a p = 0) for a set of primes of density 1 in order to conclude that f(z) = 
Sn>i °"n.q n has CM by <p. Furthermore, this implies that eip 2 = e, so ip must be a quadratic character. 

Remark. Using Ti(N), we can define the notion of a CM cusp form on all T, namely in the direct limit 
over all levels. 

The theory of ^-adic representations gives a way to generate cusp forms on Ti(N). In particular, Hecke 
characters of imaginary quadratic fields give rise to CM cusp forms. 

1.1 Cusp forms attached to Hecke characters 

To every algebraic Hecke character tp, we can attach a cusp form in a canonical way by describing the g-series 
expansion in terms of the values the character takes outside of primes dividing the conductor. More precisely, 
if ip has conductor f and co-type t on an imaginary quadratic field K, define the q-series 

/*,*(*)= E Ha)-Q NMa) (<Z = e 2 ^,Im(z)>0). (9) 

coprima S to f 

This is well-defined as a function on the ^-coordinate from the upper half plane to C. One can further check 
that it has certain symmetric properties that allow Jk,iP to be viewed as a modular form. 

Let K and ip be as above, let uj^ be the Dirichlet character of conductor Nm(f) attached to ip, defined by 

0J^(a) = ip((a)) I a* Va £ Z. 
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Note that t € Z>o due to algebraicity. Furthermore, let \k be the quadratic character attached to K, 
defined by the Kronecker symbol, 

XK(P) = (y ) , 

for odd primes p not dividing the discriminant dpc, and extended to Z in the usual way. 

Theorem 2 (Hecke |10j . Shimura |25] & [26] )■ T7ie q-series fK,ip(z) is a newform of weight t+1 and level 
d ■ Nm(f) in the space 

G S t+1 (T (d ■ Nm(f)), X K ■ w*)) C 5^.1 (T x (d • Nm(f))). 
In addition, distinct cusp forms fK,ip arise from distinct pairs (K, ip) . 

Thus, given a weight k, level N, and Nebentypus e, we can construct a basis for the space of CM cusp 
forms using the criterion in the above theorem. 

Corollary 2.1. For any positive integer r, 

fxMr ■ *) = J2 ^ ■ < ?, " Nm(a) £ S t+1 (T (N), e) <=► r ■ d ■ Nm(f) | N and X k • = e 

n 

«/ /ias oo-type t. The second equality takes place while viewing X K> kty £ is characters on Z. Thus, 
for all primes p\N, Xk(p) ■ £«ty(p) = X {p)- 

Note that the cusp forms fx,^ {r-z) have CM by xk ■ (By construction, the coefficient of q p in the g-series 
expansion of (r ■ z) is if no ideal of K has norm equal to p. Since Xk(p) = — 1 exactly when this 
holds for p, a p = Xk{p)o- p .) However, it is not at all obvious that these are all the cusp forms with complex 
multiplication in Sk(T\(N)). Using the theory of Galois representations, Ribet proved that a newform / has 
CM by an imaginary quadratic field K if and only if it arises from a Hecke character on K (see 4.4 and 4.5 

in mi 



2 The exceptional cases of n < 10 

One can try to show from the classical formulas for r n (m) that 9 n (q) are elementary when n < 12 and even. 
Instead, simply considering the dimension of the first few weighted spaces of modular forms on Ti(4) allows 
us to make the same conclusion that n = 2, 4, 6, 8, and 10 are elementary. From the Riemann-Roch formula 
specialized to ri(4), we have the following dimension formulas (see [6], [17], or [5] for proof): 

Lemma 1. The dimensions of A4fc(ri(4)) and its subspace of cusp forms for each k G Z>o are as follows: 

'0 i/fc<4 



dim c (X fe (r 1 (4))) = { HI 1 * " ™™ dim c (S k {T 1 (4))) = { k -± if k > 3 is even 



£±i if k is odd, 



^ if fc > 3 is odd. 



In particular, this implies that .Mfc(ri(4)) is entirely Eisenstein for k < 5, and iSs(ri(4)) has dimension 
1. Thus, the backwards direction of Theorem Q] follows from proving that the generator of iS>5(Fi(4)) is a 
CM cusp form. By Corollary 12. 11 the CM newform attached to the Hecke character of trivial conductor and 
oo-type 4 on Q(i) 

fowA*) = E V>(a)-? Nm(a) where </>((«)) = a 4 , 

a integral 

/ \ 

E E » 4 



4 

m>l 



\Nm(0)^m / 
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It is an element of 5s(ro(4), Xi) = $5(^1(4)), where \4 denotes the nontrivial Dirichlet character of conductor 
4. Thus, since /Qft) i/> is a newform, it is the generator of the entire cuspidal space, hence iSs(ri(4)) = 

s 5 cm (ri(4)). 

From the above argument and the fact that 9 n (q) G A4 n / 2 (ri(4)), we can therefore conclude that #2(9), 
04(q), Oeil), #8(9), and 9io{q) are elementary. 



3 n (q) is not elementary for n > 10 

The previous section illustrates that 9 n (q) is elementary if n = 2,4,6,8, and 10 using the dimension of 
■A4fe(ri(4)) for 5 = fc < 5. However, Lemma Q] also implies that £fc(ri(4)) has dimension 3 for even fc > 2 
and dimension 2 for odd k > 2. We can go further and prove the entire theorem using a dimension argument 
via the following lemma. 

Lemma 2. The dimension of iS£ m (ri(4)) is 1 if k = 1 mod 4 /or fc > 5 cmd otherwise. 

Proof. Assume iV = 4. The two Dirichlet characters on (Z/4Z) X are the trivial character 1 and the primitive 
character of conductor 4, \i- Recall that the decomposition of Sk(Ti(N)) = ® £ Sk(To(N), e) occurs over 
all characters defined modulo N such that e(— 1) = (— l) fc , thus 



5 fc (Ti(4)) 



5 fc (r (4),l) if fc is even 
5 fe (r (4),X4) if fc is odd. 



Let e denote the correct Dirichlet character such that <Sk(Fi(4)) = 5fc(r (4),e) depending on the parity of 
fc. We first want to understand the CM subspace S^ m (To(4), e) by considering the newforms which generate 
it according to Corollary 12. II 

For a newform fK,ip(r • z) G <Sji m (ro(4), s), the imaginary quadratic field K must have discriminant da 
such that lift- | 4. Thus, if = Q(i) since there is no field of discriminant —2. Note that \k is a nontrivial 
Dirichlet character of conductor 4, so \k = X4- 

Since the discriminant of K is —4, we are left with finding all Hecke characters ip on K of conductor 1 
such that w^, = 1 if the oo-type t of t/> is even, and w^, = %4 if t is odd. By definition, if a G _ftT x then since 
(a) = (ua) for all u G O^, 

= a* = {uaf = ^i((ua)). 

Since O^- = (i), we can conclude that 4 | t, since if u = ±i, i* = 1 if and only if t = mod 4. In 
particular, there are no Hecke characters of odd oo-type, hence 5^ m (r (4),e) is trivial for even fc. In 
addition, 5^ m (ro(4), e) can only be non-trivial when fc = 1 (mod 4) and fc > 1. 

We now prove that there is a unique Hecke character with the properties described above for each oo-type 
t E Z such that 4 | t. Let t be fixed such that t = (mod 4). Then xk ■ kty, = X4 f° r an y hence the 
attached character w$ must be trivial. Thus, since Ok = is a PID, the algebraic Hecke characters on 
K are exactly those defined as 

tpt ■ (") 1 — > ci Va G K x where t G Z and 4 | t. 

Using Corollary 12.11 we conclude that the dimension of S'^ m (ro(4), e) C 5fc(ri(4)) is 1 if fc = 1 mod 4 and 
fc > 5 and is otherwise. □ 



3.1 Proof of Theorem [T] via a dimension argument 

We want to show that one must calculate the coefficients of a non-CM cusp form in order to calculate the 
representation number for each integer by the quadratic form x\ + ... + x\ when n > 10 and even. To do 
this, we first construct a canonical basis for the space £„/ 2 (ri(4)) ® 5°%(Fi(4)). 



G 



Consider first if the weight § = k > 2 is even. Then Mk(^i(^)) = £fc(r (4), 1 4 ). By Lemma [TJ the 
dimension of the space is 3. 

Let I be a prime. For / e £fc(Li(4)), the associated 2-dimensional reducible ^-adic representation will 
have the form pi = e\ @£2x\~ 1 1 where the Dirichlet characters e± ■ £2 = I4 (see [21], Pp. 28), and \£ denotes 
the ^-adic cyclotomic character. Furthermore, since the product of the conductors of £\ and £2 must divide 
the desired level N = 4, both characters must be trivial. Consider the representation 



fe-i 



P = l®X e 

The L-series attached to this represenation is 



GL 2 



P P m>l \d|m / 



Taking the inverse Mellin transform gives us the q-series 



E[z) 



a 



fc-i 



(10) 



In order to calculate ao, we consider the functional equation for L-series and check that its Lambda function 
A(s,p) at s = has residue — ao = — C(l — in terms of the Riemann zeta function. Thus, ao = — Ijr 

where bk is the fcth Bernoulli number (see Section 2]). 

Using the classical converse theorem of Weil, we conclude that E(z) lies inside £fc(SL 2 (Z)) C £ k (Ti(4)) 
since the conductor of the trivial character is 1 (see [T7], §4.7, and [TS], Ch. 9). However, this implies that 
the q-series E(q 2 ) and E(q 4 ) are also in £fc(Li(4)), and furthermore, they are linearly independent. These 
are all eigenforms, and since the dimension of £fc(ri(4)) is three, B — {E(q), E(q 2 ) 7 E(q 4 )} is a basis. 

Since 9 n (q) is a modular form of weight k = n/2, we are considering the case when n = mod 4 for n > 8. 
Thus, if 9 n (q) is elementary, then there exist constants a, 6, c G C such that 

9 n {q) = a ■ E(q) + b ■ E(q 2 ) + c ■ E{q 4 ). 

Consider the determinant of the matrix of coefficients of the 4 q-series: 



det(M) = 



2n 
1 





4Q) 
1 + 2"/ 2 - 1 

1 





V3/ 






16(2) +2n 

2 n/2-l + 4 n/2-l 

1 + 2"/ 2 - 1 

1 



{coefficients of d n ) 
(coefficients of E(q)) 
(coefficients of E(q' 2 )) 
(coefficients of E(q 4 )) 



The determinant of M is zero if and only if there is a linear dependence amongst the coefficients, i.e. if the 



coefficients of q, q 2 , 
determinant gives 



q , and q of 9 n can be written in terms of coefficients of the forms in B. Solving for the 



det(M) 



2n- (1 + 3-?- 1 ) -8 



-2n-2n-3 = " 



n(n - l)(n - 2) 
6 ' 



Note that the growth of this function is dominated by —In ■ 3 a" -1 , and for n > 8, det(Af) as a function 
on n is monotonically decreasing. The determinant is when n = 4 and 8, and negative for n — 12, thus 
det(M) is nonzero for all larger n. This implies that 9 n ^ £„/2(ri(4)) for n = (mod 4) when n > 8, hence 
by Lemma [21 this implies that 9 n is not elementary for n > 8 when n = (mod 4) (as there are no CM 
cuspidal forms of even weight). 
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Now assume n = 2 (mod 4), i.e the weight of n {q) is odd. When n = 6 (mod 8), there is no CM subspace 
of <S„/ 2 (ri(4)) since k = § = 3 mod 4. We hrst check that 9 n £ n / 2 (Ti(4)) for all n = 2 mod 4 and n > 6, 
which will reduce the problem to whether there is a contribution by a CM cusp form when n = 2 mod 8. 

When n = 2 (mod 4), the £k(Ti(4)) is a two-dimensional subspace by Lemma [TJ Since the weight k is 
odd, £fc(ri(4)) = ffc(r (4), X4). Thus, we want to construct 2 linearly independent Eisenstein eigenforms of 
level 4 with Nebentypus %4- 

Let I be a prime. Similar to before, a 2-dimensional reducible £-adic representation associated to a basis 
element of £/c(Lo(4), xa) will have the form e\ © £2X£~ 1 > where the characters £1-62 = X4 an d the product of 
their conductors divide the level (see [21] . Pp. 28). Thus, we have the following two possible representations: 

— ► GL 2 (Qi), 



Pi 
92 ■■■ 



X4 



1 . 



GL 2 



The L-series attached to these representations are 



L(s,pi) 



L(s,p 2 ) 



n 1 - X4 (p>- s n 1 ^-ip-- n ( n x4 ( d ) dfc 1 m . > 

p p m >l \d|m / 

tt — ^ — tt — — = y (y X 4(d)d fe - 1 1 — . 

11 1 _ X4 (p)pfc-lp-« I ^ W J 771 s 



Applying inverse Mellin transforms to the above L-functions produces two g-series 



Ex(q) 



ao,i 




ik-i 



(11) 



E 2 {q) 



(12) 



One can calculate ao,i and ao,2 by considering the functional equations for the L-series and calculating the 
residues of the associated Lambda functions A(s,p±) and A(s, p 2 ) at s = 0. One can check that —00,1 = 
— L(— 2,X4)/2, in terms of the L-series associated to the Dirichlet character X4, and — ao,2 = 0. Thus, 



b £ 4 

ao,i = — and a , 2 = 0, 



where 



t 



t ■ e 



•M 



1 



1 



OO 

Vb»- 



fc=0 



is the fcth Bernoulli number associated to the nontrivial Dirichlet character of conductor 4 defined by the 
above equation (see [35], Ch. 5). 

The converse theorem of Weil implies that both E\{q) and E 2 (q) are distinct eigenforms on Li(4) of weight 
k, hence the two series are linearly independent (see [T7], §4.7, and [15] . Ch. 9). By construction, Ei(q) and 
E 2 {q) £ £fc(ri(4)), thus we take B = {Ei{q) , E 2 {q)} as a basis for the Eisenstein space, £ n / 2 (Tx(4)). 

To show 8 n ^ £ n / 2 (Ti(4)), we check if the following matrix of coefficients for q, q 2 , and q 3 has non-zero 
determinant: 



2n 
1 
1 



4(2) 
2«/ 2 — 1 

1 



-1 



det(M') 

Solving for the determinant, 

det(M') = 2n h^ 2 ' 1 + l) (l - T /2 -^ - 8 

= (-4n 2 + 6n)+2"/ 2 - 1 -(-^n 3 + 



V37 

g»/2- 



1 - 3"/ 2 - 



{coefficients of On) 
{coefficients of E\{q)) 
{coefficients of £2(9)) 



r») + 3"/ 2 - 1 • (4n 2 - 6n) + 6"/ 2 - 1 • (-2n) . 



4n 2 - 
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Note that for n > 6, det(M') as a function on n is governed by the growth of — 2n ■ 6™/ 2_1 . One can check 
that for n = 6, det(M') = 0, and for n > 6, det(M') < and is monotonically decreasing. We conclude that 
for n = 2 mod 4 and n > 6, 9 n is not a linear combination of Eisenstein series. By Lemma [21 this implies 
that for all positive n = 6 mod 8 except n = 6, 9 n is not elementary. 



It remains to show that n is not elementary when n = 2 mod 8 and n > 10. Since the weight k = § 
1 mod 4, the space iS^(Ti(4)) is generated by a single CM cusp form C(q) constructed below. 

Let ip denote the algebraic Hecke character unramified away from 2 on Q(i) with oo-type j — 1. Viewing 



it adelically, ip acts on C x by sending zi-yz 1 fc , and on primes p £ Ok, "4> sends p H> p k 1 , where k = ^ as 
usual. The i-series attached to this character is 



(1 



+ (i + 



k-\< 



p=3(4) 



1 



n 



p = l(4) (1 
P=PP 



S )(l 



-fc-1 



, oez[i' 
\Nm(c)) = 



Applying an inverse Mellin transform gives the g-series 

C(g) 



m>l 



\Nm(0)=m / 

g + (-4) fe -Y + 2 fe - 1 ( Z 4 + 



(13) 



By Weil's converse theorem, this is a cusp form on Li(4) of weight k, and by construction, it has complex 
multiplication by the Dirichlet character \i attached to Q(i) (see [17], §4.7, and [15], Ch. 9). Thus C(q) € 
iS£ m (ri(4)) and in fact, generates the space. 

Thus, 8 n is not elementary if the determinant of the matrix of coefficients for q, g 2 , g 3 , and g 4 , 



det(M") 



has non-zero determinant. 

A straightforward calculation demonstrates the fact that det(Af') as a function on n has negative growth 
on the order of C(12™/ 2 - 1 ) for n large enough. For n = 10, det(Af') = 0, but for n = 18, det(M") = 
—439, 038, 812, 160. It is therefore easy to check that det(M") is nonzero for all n > 10 such that n = 2 mod 8. 
This proves the theorem using the dimension arguments from Lemmas [1] and [2] □ 



2n 


4© 


8Q) 16f 4 l ) +2n 


(coefficients 


of8 n ) 


1 


2^/2-1 


_l + 3«/2-i 4"/2-l 


(coefficients 


o}E x (q)) 


1 


1 


1 - 3"/ 2 - 1 1 


(coefficients 


of Mi)) 


1 


(_4)V 


2^ 


(coefficients 


ofC(q)) 



4 n (q) is not elementary for n > 10, redux 

Using the construction of a basis for the CM cuspidal subspace in Corollary 12. 11 we can show that the Hecke 
action of the cuspidal part of 9 n (q) does not coincide with the action on elements of 5n m (F 1 (4)) at the prime 
p = 3 for n > 10. This allows for another proof of Theorem [1] We first describe Hecke operators which 
vanish on CM cusp forms. 
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Lemma 3. For any f £ <S£ m (Ti(4)) with k > 1, the Hecke operator T p (f) = for any prime p = 3 mod 4. 

Proof. Recall that Corollary 12.11 implies that S£ m (ri(4)) is generated over C by forms fq(i)^{z) defined by 
©, where ip is a Hecke character with co-type t divisible the level N = 4. In order for any fqu) t ^(z) to have 
CM by the quadratic character xa, the pth coefficient must vanish, i.e. a p = for all inert primes p, for all CM 
eigenforms. Thus, for any prime p = 3 (mod 4), T p (f) must necessarily vanish for any / £ S" n (Ti(4:)). □ 

Remark. Note that this result can be easily generalized to any level N: for any / £ S^, m (Ti(N)) with 
k > 1, T p (f) = for any prime which is inert in all imaginary quadratic fields K with discriminant — d 
where d\ N. It is also possible that for a larger set of primes, T p (f) vanishes (e.g. for Ti(144), T p (f) = 
for all primes p = — 1 mod 12 or equivalently, all primes which are inert in only Q(i) and Q(\/— 3), see [24], 
Lemma 1). Furthermore, the converse to this statement is also true: any cusp form / £ <Sfc(ri(4)) such that 
T p (f) = for aUp = 3 mod 4 is an element of S£ m (ri(4)) (see [23], Corollary 15.2). 

Thus, if we show that the part of 9 n which vanishes on all three cusps of T\ (4) has a nonzero coefficient for 
q p for some p = 3 mod 4, then 9 n is not elementary. Since 9 n e A / J„/ 2 (ri(4)) = £ n / 2 (ri(4)) © 5 n / 2 (ri(4)), 
we can decompose 9 n as the sum of an Eisenstein series and cusp form, i.e. write 9 n {z) — E n /2(z) + f n /2(z)- 

We use a result derived from the Siegel-Weil formula describing the coefficients of E n / 2 (z) (see [30], pp. 
373-376). Note that because we are focused on Hecke operators of odd primes, we can focus on the odd 
coefficients of the Fourier expansion. In order to write down formulas for the odd coefficients, we must first 
understand certain constants. 

Let et denote the fcth Eulcr number and bj denote the jth Bernoulli number, which are defined by the 
following identities 



The magnitudes of these numbers are related to values of L-series associated to the primitive Dirichlct 
characters of conductor 1 and 4: 

2 • j 1 2 2fe+3 • fc 1 

b J = To ii ' if i > even > and e k = ln TuTT ■ L ( k + 1, Xa) if fc > even. 

Theorem 3 (Shimura [55] and 28 ). If E n / 2 (z) — 1 + Em=i f™?™ ^ ^n/2(ri(4)) denotes the q-expansion 
of the Eisenstein part of 9 n /i[q), then for odd m. 



This result comes from using local computations on the genus and the mass of the standard lattice Z ra to 
understand the "dominant" growth factor c m of r n (m). To demonstrate that r n (m) is generally governed 
by the growth of c m , one shows that c m — r n (m) are coefficients of a cusp form, thus 1 + 53m=i c m<f n is the 
Eisenstein part of the 9 n (q). 

Lemma 4. If f n /2{ z ) denotes the cuspidal part of 9 n {z), then the 3rd coefficient as{f n /2) 7^ if n > 10. 

Proof. One can easily compute that r„(3) = 8(") = 4n ("~i)("~ 2 ) _ We can then calculate the third coefficient 
of the g-series expansion of f n /2{ z ) = 9 n (z) — E n /2(z) using this equality: 






if n = (mod 4) 



00 



00 00 



fn/2(z) = ^ am< i 



rn(m)q m - ^ c m q m G 5 n/2 (r 1 (4)). 



m— 1 



m— m— 
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Thus, taking the difference of c 3 and r„(3) gives the coefficients of f n /2{ z ) 
' An(n - l)(n - 2) 



a 3 



4n(n — l)(n — 2) / // 



+ (raH 2 "'" 1 - 14 ©H 1 - 3 "' 2 ") ,fns2 (mod4) 

1+3 n/2-l> 



1/2 1 



1 - 2"/2 



if n = (mod 4) 



Note that the values of |bj| and |ek| for even j and k have bounds given below (see [13]), 



2 \2-Ke J 

Thus, we can bound a 3 as a function on n by 



and y <|ek| - 



a 3 (n) > 



4ra(7i- l)(n- 2) ( 1 2w ( ire 



2 V n/2-1 

■1/2 



n/2-l N 



• ( 2™/ 2 " 1 — — 3 n/ 



a 3 



4n(n-l)(n-2) /n /T" /47re\ ™ /2 \ /l + 3 n / 2 - x \ 



It is easy to check that for n large enough, the non-polynomial part of the lower bounds has negative 
growth on the order of 0{n~ n / 2 ) which is dominated by the positive growth of the polynomial portion. 
Some asymptotic calculations further demonstrate that these functions on n are monotonically increasing 
for n > 18 and 20, so in particular 0,3 is nonzero for n > 10 as seen in the following table.) 



n 


c 3 


r»(3) 


a 3 = r n (3) c 3 


4 


32 


32 





6 


160 


160 





8 


448 


448 





10 


960 


960 





12 


1952 


1760 


-192 


14 


189280/61 


2912 


-11648/61 


16 


70016/17 


4480 


6144/17 


18 


1338240/277 


6528 


470016/277 


20 


157472/31 


9128 


125248/31 



Thus, the cuspidal part of 6 n has a nonzero coefficient for q 3 in the Fourier expansion. The first g-series 



coefficient of T 3 (s n / 2 ) is ^3 which is nonzero for n > 10. 



□ 



Combining the two lemmas, we conclude that the cuspidal part of 9 n does not lie in the CM subspace 
of iS„/ 2 (ri(4)) for n > 10, hence 8 n cannot be elementary for these cases. Theorem [1] again follows from 
showing that 8 n is elementary for n = 2, 4, 6, 8, 10 from Section [2j 



5 Motivation for definition of elementary modular forms 

While the notion of an elementary modular form coincides with the property of having lacunary cuspidal part 
as well as being attached to potentially abelian ^-adic representations of Gal(Q/Q), the motivation behind 
the definition of "elementary" is computational: an elementary modular form should have Fourier coefficients 
which have formulas that can be computed in polynomial time in log(m) (for the coefficient of q m ). Since 
the mth coefficient of an Eisenstein series is usually given as a sum over the positive divisors of m, one can 
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weaken this by assuming factorization of m is given when calculating the coefficient of q m in log-polynomial 
time. For 0„, this occurs for even n < 12, demonstrated by the formulas given in these exceptional cases, 
and these allow for efficiently computing r„(m), even for 0io which is not purely Eisenstein. To demonstrate 
that modular forms on Ti(4) that involve a non-CM cusp form are not "elementary," we consider the case 
of n = 12. 

5.1 n = 12 

Note that the subspace of A^6(ri(4)) consisting of Eisenstein series has dimension 3, and the subspace of 
cusp forms has dimension 1 by Lemma [TJ Furthermore, by Lemma [3J <Sg'"(ri(4)) is trivial, and we can in 
fact use a well-known cusp form for the basis element of iSe(ri(4)): 

oo 

>/A(2z)= ?7 12 (2z) := q J| (1 - q 2n ) 12 = q [l - I2q 2 + 54g 4 - 88q 6 - 99g 8 + ...] 

n=l 

The Dedekind eta function 77(2) is a cusp form on Ti(2) of weight 1/2, hence ^/A(2z) 6 A^6(Fi(4)). It is a 
cusp form without complex multiplication. 

From above and (JTUJ), take {E(q), E(q 2 ), E(q 4 ), VK(q 2 )} with k = 6 as a basis for A^ 6 (Li(4)). We then 
calculate the constants a, 6, c, d € C such that 

6>i 2 (g) = a ■ E{q) + b ■ E(q 2 ) + c ■ E(q 4 ) + d ■ VA(q 2 ). 

Note that by Theorem [TJ d must be nonzero. The hrst few coefficients of the g-series in the basis are as 
follows: 

E(q) = -^+q + 33-q 2 + 2U-q 3 + l057-q A + ... 
E(q 2 ) = -^+q 2 + 33-q i + ... 

VA{q 2 ) = q-12-q 3 + ... 
Since we know that the first few coefficients ri2(m) are 

0i 2 (g) = 1 + 24 ■ q + 264 • q 2 + 1760 • g 3 + 7944 ■ g 4 + 

thus a — 8, b = 0, and c = —512, d = 16. 

Here, we can conclude that it is necessary to calculate the coefficients of rj 12 in order to calculate 0i2(n). 
Furthermore, the converse is true. 

From Serre's point of view, 77 12 is not lacunary, i.e. it is not a CM form, thus the coefficients cannot be 
given via characters of a CM field (see [23], §15). Furthermore, it has a positive density of non-zero coefficients 
in its g-series that contribute to 0i2 and are hard to compute (see [24 ). The formulas for n = 2,4,6,8, and 
10 illustrate that r„(m) can be calculated efficiently if the prime factorization of m is known (thus for a 
prime p, r n (p) can be computed in log(p)-time). In the case of n = 12, this is not enough information as 
there is no analogous description of the coefficients of rj 12 (2z) in terms of divisors. One must know how to 
compute the Ramanujan t function in order to calculate ri 2 (m). 

5.2 Making non-elementary 9 n (q) computable 

Recently, Bas Edixhoven, Jean-Marc Couveignes, and Robin de Jong have proven that if / = ^ a n q n is a 
modular form on SL2(Z), then even when there is not an elementary formula for the Fourier coefficients, 
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a p for p prime can be computed in time polynomial in the weight k and log(p) assuming GRH (see [5]). 
This implies that the prime Fourier coefficients r(p) of A can be calculated in polynomial time with respect 
to log(p). Furthermore, Peter Bruin gives a description of a probabilistic algorithm in time polynomial 
to k and log(p) which under the assumption of GRH, computes the Fourier coefficient of q p of eigenforms 
of level 2N' where N' is squarefree (see [5]). This includes calculating r n (m) for all even n as discussed 
here, whether or not 9 n has an elementary formula. In particular, since 9 n is not elementary for all even 
n > 10, and therefore does not have a nice formula, such an algorithm is clearly necessary. Note that 
for elementary 9 n , the formulas given here compute coefficients faster than Bruin's algorithm. Currently, 
Bruin's algorithm is largely theoretical and future work will be done to make it more practical. Other than 
understanding classical arithmetic functions, these algorithmic results are useful in computing eigenvalues 
of Hecke operators or equivalently, coefficients of eigenforms (see [7J, [3], and [5]). 
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